Phonon spectrum in a nanoparticle mechanically coupled to a substrate 
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We calculate the vibrational density-of-states in an insulating nanoparticle that is in weak me- 
chanical contact with a semi-infinite substrate. The work is motivated by a recent experiment by 
Yang et al, where the low-energy phonon density-of-states of Y2O3 nanoparticles doped with Eu^^ 
was measured. Preliminary results presented here, based on the conventional quasiparticle-pole 
approximation for the phonon propagator, are in reasonable agreement with experiment. 



I. INTRODUCTION 

As is well known, the vibrational spectrum of an iso- 
lated nanometer-scale crystal is discrete. For a spheri- 
cal nanoparticle of diameter d and bulk transverse sound 
velocity Vt, the lowest internal vibrational mode, called 
the Lamb mode, has a frequency of about 2TTVt/d. Be- 
low this frequency no internal vibrational modes exist, 
and the phonon density-of-states (DOS) vanishes. Any 
property of the nanoparticle that depends on this vibra- 
tional spectrum, such as its thermodynamic properties 
or electron-phonon dynamics, will therefore be very dif- 
ferent, especially at low energies, when compared to bulk 
crystals made of the same material. 

In an interesting recent experiment, Yang et al. |^] 
measured the DOS of Y2O3 nanoparticles with a size dis- 
tribution ranging from 7 to 23 nm in diameter. Nanopar- 
ticles of size, say, 15 nm, cannot support phonons with 
energies below about 10 cm~^. The phonon DOS was 
obtained by measuring the nonradiative lifetime of an ex- 
cited electronic state of Eu'^+, and at 3 cm"-"^ was found 
to be about 100 times smaller than that of a bulk Y2O3 
crystal. 

In this paper we propose and investigate a mechanism 
to explain the observed low-energy DOS. Several phonon- 
broadening mechanisms could be responsible for this ef- 
fect. For example, anharmonicity causes the modes to 
broaden, leading to a DOS at low energy. However, an- 
harmonicity is ineffective at low energy and an estimate 
of the anharmonic broadening showed that the resulting 
DOS is much smaller than that observed ||^. Adsorbed 
molecules or "dirt" on the outside of the nanoparticle 
could also be involved. A more likely broadening mech- 
anism results from the fact that these nanoparticles are 
not isolated, but rather are in contact with each other 
or some support structure. This contact enables the 
nanoparticles to couple to an environment with a con- 
tinuous spectrum at low energy. We believe that it is 
this physical coupling to the environment that broadens 
the vibrational modes enough to explain the observed 
DOS. 

Our preliminary results, presented here, are favor- 
able. However, the calculation has only been done within 
the conventional quasiparticle-pole approximation for the 
nanoparticle phonon propagator which, to leading or- 
der in perturbation theory, is equivalent to using Fermi's 



golden rule to obtain vibrational-modc lifetimes and as- 
suming Lorentzian line shapes. A more accurate calcula- 
tion of the nanoparticle DOS, which requires the solution 
of the Dyson equation, will be presented in future work. 



II. MODEL AND METHOD 

Our model consists of a single isotropic elastic sphere 
of diameter 10 nm, representing the nanoparticle, con- 
nected to a semi-infinite isotropic elastic continuum lying 
in the xy plane and extending to infinity in the negative 
z direction. For simplicity, we take the substrate and 
the nanoparticle to be made of the same material. We 
model the contact between the two by a weak harmonic 
spring, corresponding to the situation where the nanopar- 
ticle and substrate are connected by only a few atomic 
bonds or by a small "neck" of material. The Hamiltonian 
for the system is 

H — TiLUjo^jaj + TiLOib\bi 
J I 
1 2 

+ 2 -^['"^(''o) part - Mz(ro)sub] , (1) 

where the a,j and a\ are annihilation and creation op- 
erators for phonons in the nanoparticle, and the index 
J runs over all the modes of the nanoparticle. The hi, 
b\, and / correspond to the substrate phonons [Q. The 
spring constant K is taken to be of the order of an atomic 
bond strength of the material. Uz(ro)part and Uz(ro)sub 
are the z components of the phonon displacement field 
of the nanoparticle and the substrate respectively, eval- 
uated at the point of connection Tq. 

The displacement field for the nanoparticle can be ex- 
panded as 



u(r, t)part — 




[aj*j(r)e-*"-'*+H.c.], (2) 



where the ^'./(r) are the nanoparticle vibrational eigen- 
modes, normalized according to 

/ d\<b:j{v)-<i>j,{v) = 5j.r, (3) 
Jv 



1 



and p is the mass density. V is the volume of the nanopar- 
ticle. The nanoparticle's vibrational eigenmodes were 
found by using a method similar to the one used by Lamb 

i- 

Similarly, the substrate displacement field is given by 



u(r, t) 



sub 



[6/f/(r)e 



— iujjt 



H.c], (4) 



with the f/ normalized as above. The vibrational eigen- 
modes of the substrate were calculated following Ezawa 
§. 

To determine the vibrational DOS we calculate the re- 
tarded Green's function for the nanoparticle, 



(5) 



Evaluating this perturbatively to second order in the 
spring constant K , and using the quasiparticle-pole ap- 
proximation 1^ (which replaces the energy-dependent 
phonon self-energy for each mode J with it's on-shell 
value), we find the energy damping rate for mode J to 
be given by 



hp ujj 



*}(ro)| 



(6) 



where 



and mass density p = 2.3 g cm~'^. For a spring stiffness K 
equal to 1.0 X 10^ erg cm~^, we obtain the phonon DOS 
shown in Fig. |l|. 
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FIG. 1. Vibrational DOS in a 10 nm Si nanoparticle, 
weakly coupled to a substrate. The DOS rapidly ap- 
proaches the bulk Debye law. 



N,{lu) = ImDr(ro,ro,u;) 



(7) 



is the phonon spectral density at the surface of the sub- 
strate and Dl^{r,r'uj) is the transform of the substrate 
propagator. Eq. can also be derived from Fermi's 
golden rule. The quasiparticle-pole approximation pre- 
dicts Lorentzian-broadened line shapes. This leads to a 
nanoparticle DOS given by 



5(e) 



(e — hujj)^ + ^ '''j 



(8) 



The usual thermodynamic DOS (number of states per 
unit energy per unit volume) is given by g(e)/V, where 
V is again the volume of the nanoparticle. 



III. RESULTS AND CONCLUSIONS 

For simplicity we assume the nanoparticle and sub- 
strate to be made of Si, thus enabling us to obtain the 
surface spectral density from Appendix B of Ref . . 
We treat Si as an isotropic elastic continuum with longi- 
tudinal and transverse sound velocities 



vi = 8.5 X 10^ cm s \ 
vt = 5.9 X 10^ cm s"\ 



Our calculation includes all modes below a cutoff fre- 
quency Wmax of 100 cm^^. The low-energy DOS depends, 
to some extent, on our choice of Wmax- This sensitivity 
is an artifact of the quasiparticle-pole approximation, or 
equivalently, a consequence of assuming Lorentzian line 
shapes. Because at this stage we are only interested in 
whether our result agrees with the experiment of Ref. 0] 
at the order-of-magnitude level, we will not discuss the 
weak Wniax dependence any further. 
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FIG. 2. Low-energy DOS. 
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Fig. H shows the low-energy DOS up to about 15 
cm~^. The large peak on the right is the five- fold- 
degenerate Lamb mode. At 3 cni^^ the DOS is found 
to be 5.8 X 10~^ states per cm^^, or about 3.2 x 10^^ 
times the bulk DOS il 



E 



V = 2.0x 10--^ E-" 



states 

_ 1 ; 



(10) 



where E is the energy in cm^^. Experimentally, the ratio 
of nanoparticle to bulk DOS was found to be approxi- 
mately 7.4 X 10^^, about 20 times larger than our result. 

This is a good result considering the simplicity of our 
model. However, the quasiparticle-pole approximation is 
probably inaccurate much below 10 cm^^, where the de- 
viation from Lorentzian line shapes becomes important. 
However, the order-of-magnitude agreement does suggest 
that we have correctly identified the relevant broadening 
mechanism in these nanoparticles. An accurate calcula- 
tion of the nanoparticle DOS, based on the solution of 
the Dyson equation, is in progress. 
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